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ABSTRACT 

A gyro-based filter variation on the standard QUEST attitude determination 
algorithm is applied to the Cosmic Background Explorer (COBE). Filter 
QUEST is found to be three times as fast as the batch estimator and slightly 
more accurate than regular QUEST. Perhaps more important than its speed 
or accuracy is the fact that Filter QUEST can provide real-time attitude solu- 
tions when regular QUEST cannot due to lack of observability. Filter QUEST 
is also easy to use and adjust for the proper memory length. Suitable applica- 
tions for Filter QUEST include coarse and real-time attitude determination. 


* This work was supported by the National Aeronautics and Space Administration (NASA)/Goddard Space 
Flight Center (GSFC), Greenbelt, Maryland, Contract NAS 5-31500. 
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1. INTRODUCTION 


Filter QUEST is a new version of the familiar QUEST attitude determination software. Be- 
cause it makes use of gyro data, it promises to combine the speed of QUEST with the accuracy 
of gyro-based estimators. Just how fast and how accurate is seen from experience with 
Cosmic Background Explorer (COBE) data. 

2. BACKGROUND ON QUEST 

QUEST is Shuster’s implementation of the Davenport q-algorithm for attitude determina- 
tion (Reference 1). Since its introduction in 1978, it has been a standard part of the attitude 
ground support systems built for three-axis spacecraft at the NASA Goddard Space Flight 
Center Flight Dynamics Facility. Explanations of the algorithm and its implementation are 
available elsewhere (References 2, 3, and 4), but a few remarks about QUEST that are rele- 
vant to the following discussion are included for convenience. 

First, QUEST uses only unit vector observations expressed in body coordinates. These are 
denoted here by w ( , where i indicates the particular observation. The corresponding refer- 
ence vectors in inertial coordinates are Vj . Each observation is also given a scalar weight aj, 
which is normalized so that the sum over all observations equals one. 

N 

I a, = 1 (1) 

i = 1 

Second, an intermediate quantity called the attitude profile, or “B-matrix,” is computed as 
follows: 


B =* ^ aj Vj wF 

i = i 


Several secondary quantities are also defined: 


z 


B23 - B32 
B31 - B13 
B12 - B21 


( 2 ) 


(3) 


S = B + B t 


(4) 
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a = trace[B] 


( 5 ) 


Third, the attitude quaternion q is then the solution to the eigenvalue equation 

Kq = Aq ( 6 ) 

where the matrix is constructed from those secondary quantities: 

•r _ (S - Ol) Z I (7) 


3. FILTER QUEST 

Filter QUEST is a natural generalization of the QUEST attitude determination algorithm 
(References 5 and 6). QUEST is limited to using “single-frame” observations from the same 
time or a constant attitude. Filter QUEST, however, can accept observations from different 
times or attitudes so long as it has the incremental rotation angle between those times. This 
information typically comes from gyros. In the past, systems have been built that propagate 
the unit vector observations to a common time before providing them to QUEST for “batch” 
processing. Filter QUEST differs in that it handles that preprocessing internally. Filter 
QUEST also differs in that it is a sequential estimator with a simple “fading” memory rather 
than a batch estimator. In Filter QUEST, the B-matrix is corrected for spacecraft rotation 
between observation times tj and tj.i< Here Bj/j denotes the value of the B-matrix at time tj 
based on observations up to and including those at time tj: 

Bj/j_i = aO (tj, tj_i) Bj_i/j_j (^) 

The memory length parameter a is a number between zero and one. a equal to one implies 
infinite memory, and a. equal to zero produces the standard QUEST algorithm. $ is the 

attitude propagation matrix computed from the angular increment vector 9 for the time in- 
terval (tj, tj.i). It is given by the expression 


<I> (tjtj_i) = cos 9 I + (1 - cos0) 99 T - sin 9 [0 X ] 


( 9 ) 


where 9 is the magnitude and 9 


is the direction of 9 


9 = 99 
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and 




-#3 0 2 

0 -e 

0 , o 


At each observation time, the B-matrix is updated using the new pair of vectors v ( and w, 


B ./i = Bi/i-i + aj Vj wX 


( 10 ) 


Attitude is computed at selected times using the standard QUEST algorithm. 

There are three reasons for interest in Filter QUEST. First, because it uses gyro data, it can 
provide more accurate real-time attitude solutions than regular QUEST. It can also extend 
observability into periods such as orbit night, when single-frame methods may be useless. 
Second, because QUEST is fast, it is hoped that Filter QUEST will provide accuracies ap- 
proaching those of the batch estimator with less CPU time. This is useful for producing 
24-hour attitude histories, as in definitive processing. Third, the simplicity and familiarity of 
Filter QUEST make it a suitable “training” filter. Timing is accomplished by adjusting the 
single parameter a , which can always be set equal to zero to return to the original QUEST. 
Because QUEST is so well known, any improvement is also likely to be accepted. 

To evaluate its performance with real data. Filter QUEST was built into the Coarse Attitude 
Determination Subsystem (CADS) of the Cosmic Background Explorer (COBE) Flight 
Dynamics Support System (FDSS), where it replaced regular QUEST. Tests were then made 
to see how much attitude error could be eliminated by increasing a and to determine an opti- 
mal value for a that would minimize error while avoiding long-term divergence. For these 
tests, the Fine Attitude Determination Subsystem (FADS) solution computed using Diffuse 
Infrared Background Experiment (DIRBE) star observations was used as a reference. FADS 
is a batch estimator that uses gyro data to provide the most accurate COBE attitude solutions. 
The Filter QUEST and these DIRBE batch solutions were then compared in the Quality As- 
surance (QA) subsystem, where the plots and statistics on their differences were computed. 
Finally, timing comparisons were made for the batch, regular QUEST, and Filter QUEST 
estimators. 


4. THE EFFECT OF MEMORY LENGTH 

The length of the fading memory is controlled by the parameter a , which determines the 
effect of the current observation on the current attitude estimate. Although a is constant, the 
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■'effective” weight, which is the ratio of the current observation weight to the total weight, is 
not. If all observation types have the same weight d , the sum of these weights at time tj is 


i-l 

W, = a X « J 

j=o 


a(i - «') 
l - a 


(ID 


The ratio of the current weight to the total weight is 


1 > 


a 

wi 


1 - q 
(1 - a) 


1 -a 


( 12 ) 


This number decreases from one and converges to (l - a) . Thus, the fading memory filter is 

more responsive to early observations than later ones. This is particularly so a is set high and 
there is no a priori solution. In this case, the solution can lock up quickly based on a tew 
inaccurate measurements and then take a very long time to converge to the correct attitude. 

The error in a Filter QUEST solution is taken to be the difference between that solution and a 
DIRBE batch solution. For the short timespan shown in Figure 1, the error is random with 
amplitude consistent with the 0.5 degree Sun sensor resolution. The size of the error corre- 
sponds to the spread of the difference curve or the root of the mean squared (rms) difference. 
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Figure 1 . Filter QUEST Solution Error 
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As the memory length parameter a increases from zero to one, the rms difference decreases, 
as shown in Figure 2. The rms difference does not go to zero because the observation errors 
are not completely random. 



Memory Length Parameter 


Figure 2. Error Standard Deviation as a Function of Memory Length 


It would seem then that the memory length should be made as large as possible. The adverse 
effect of doing so becomes apparent over longer timespans. As shown in Figure 3. the differ- 
ence can increase with time and exceed the one degree (3a) COBE accuracy goal. Figure 4 

gives a plot of the rms differences for full-orbit (100-minute) solutions as a function of a . 
There is no sharply defined value of a that minimizes the error. Any value between 0.90 and 
0.99 is suitable. The determination of this value is sometimes called “tuning.” In practice, the 
true attitude is not apt to be known, and tuning must be done by examining how well the ob- 
servations are fit by the attitude solution. QUEST and Filter QUEST provide such a figure of 
merit in the output “RESIDU” g , which equals one minus the “overlap eigenvalue” X (Ref- 
erence 7). When all the observations fit the attitude solution. RESIDU equals zero. 

6 = 1 - X (13) 

RESIDU is plotted in Figure 5 for the solution of Figure 3. Just as the growing difference 
curve indicates that a is set too high, the growth of RESIDU says the same thing without the 
need for an absolute attitude reference. Figure 6 gives a plot of the maximum RESIDU 
values for full-orbit solutions as a function of a corresponding to the rms differences pro- 
vided by Figure 4. Again, there is no clearly optimal value for a . Any value between 0.90 and 
0.99 will do. Alternatively, a can be chosen such that RESIDU is as large as possible but does 
not grow over time. 


144 



1 . 50001 


v> 

0) 

G) 

k. 

O) 

<D 

T 3 


O 

w 

LU 


(/) 

LU 

o 


UL 


1 . 4000 
1 . 3000 
1 . 30001 
I . 1000 
l . 0000 
0,9000 
0 . 800 0 
0.7000 
0.6000 
0 . 5000 
0.4000 


a= 1.0 


0 . 30001— 
09.30. 


00 1 0 . 0 3 . 3Q 10.36.40 i l . 10.00 

00 .46.40 10.30.00 l 0 . 5 3 . 30 11.36.40 


Time (HH.MM.SS) 

Figure 3. Divergence of Filter QUEST Solution Over a Full Orbit 


ORIGINAL PAGE IS 
OF POOR QUALITY 


145 


V) 

o 

o 

V. 

cn 

■D 


o 

w 

w 

LU 


01 

I— 

TO 

3 

a* 

(/> 

c 

TO 

<u 

2 


o 

o 

ac 



Memory Length Parameter 


Figure 4. Root Mean Square Error as a Function of Memory Length 
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Figure 5. RESIDU as an Indicator of Divergence 
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Figure 6. RESIDU as a Function of Memory Length 


5. SPEED AND ACCURACY COMPARISONS 

Since the batch estimator and QUEST represent the current COBE capabilities, they are 
used as standards against which Filter QUEST is measured. For comparison, all three estima- 
tors are applied to a full orbit (100 minutes) of data. The memory length parameter a is set 
to 0.98 to keep the attitude from diverging over this timespan. As before, the reference is a 
batch solution computed using DIRBE star data. The batch error plotted here is the differ- 
ence between this reference and a batch solution obtained from Sun and Earth sensor data 
alone. 

The rms errors of the three solutions are shown in Figure 7. For this comparison, gyro biases 
were first solved for in the batch estimator, and the corrections were made in the Data 
Adjuster (DA) subsystem. This improved the Filter QUEST solution accuracy as well as that 
of the batch estimator. In spite of having gyro data, Filter QUEST had an rms error only 
14 percent smaller than that of regular QUEST, and still far greater than that of the batch 
solution. The reason for this small improvement is that like any sequential estimator. Filter 
QUEST cannot look ahead of the current solution time. While the batch estimator effectively 
has infinite memory and thus can average out orbital and higher frequency errors. Filter 
QUEST cannot set a to one without diverging. Unless observation errors are really random, 
or are much smaller than the propagation error, no filter compares favorably with a batch 
estimator in accuracy. Unfortunately, all nonrandom errors seem not to have been elimi- 
nated from the COBE observations. 

The speed of each estimator in producing these solutions is plotted in Figure 8. These num- 
bers are for the Flight Dynamics Facility Hitachi NAS 8063 computers and come from 
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running only the estimators. The conversion of telemetry, correction for known calibration 
errors, and generation reference data are not included. These speeds are for running interac- 
tively rather than in background mode. Background execution can be much faster, but the 
relative speeds should be similar. As expected, regular QUEST is the fastest. Addition of 
propagation to Filter QUEST cuts its speed in half. It is still, however, more than three times 
as fast as the batch estimator, the speed of which has been greatly increased since it was origi- 
nally coded. 


6. CONCLUSIONS 

In conclusion, this study suggests that Filter QUEST pays in speed for what can be a small 
increase in accuracy. For this reason, it seems more suited for use as an upgrade from regular 
QUEST than as a replacement for the batch estimator. For spacecraft with only limited peri- 
ods of complete attitude observability, however, this can mean the difference between know- 
ing and not knowing the attitude in real time. Because the Filter QUEST algorithm reduces 
to regular QUEST when the memory length parameter a is zero, there is no risk in substitut- 
ing it in all current and future applications of QUEST. 
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